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Spin polarized states in dense neutron matter with recently developed Skyrme effective interaction 
(BSk20 parametrization) are considered in the magnetic fields H up to 10^° G at finite temperature. 
In a strong magnetic field, the total pressure in neutron matter is anisotropic, and the difference 
between the pressures parallel and perpendicular to the field direction becomes significant at _ff > 
Hth ^ 10^* G. The longitudinal pressure decreases with the magnetic field and vanishes in the 
critical field 10^® < He < 10^^ G, resulting in the longitudinal instability of neutron matter. With 
increasing the temperature, the threshold Hth and critical He magnetic fields also increase. The 
appearance of the longitudinal instability prevents the formation of a fully spin polarized state in 
neutron matter and only the states with moderate spin polarization are accessible. The anisotropic 
equation of state is determined at densities and temperatures relevant for the interiors of magnetars. 
The entropy of strongly magnetized neutron matter turns out to be larger than the entropy of the 
nonpolarized matter. This is caused by some specific details in the dependence of the entropy on 
the effective masses of neutrons with spin up and spin down in a polarized state. 
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I. INTRODUCTION 

Magnetars are strongly magnetized neutron stars [l| 
with emissions powered by the dissipation of magnetic en- 
ergy. According to one of the conjectures, magnetars can 
be the source of the extremely powerful short-duration 
7-ray bursts The magnetic field strength at the 

surface of a magnetar is of about IQ^'^-IQ^-' G i,^. Such 
huge magnetic fields can be inferred from observations of 
magnetar periods and spin-down rates, or from hydrogen 
spectral lines. In the interior of a magnetar the mag- 
netic field strength may be even larger, reaching values 
of about 10^^ G @, @. Under such circumstances, the 
issue of interest is the behavior of neutron star matter in 
a strong magnetic field 

A realistic description of neutron star matter should 
include, at least, neutrons, protons, electrons and muons 
subject to the charge neutrality and beta-equilibrium 
conditions. The magnetic field then influences the sys- 
tem properties through Pauli paramagnetism as well as 
via Landau quantization of the energy levels of charged 
particles. Nevertheless, because the neutron fraction is 
usually considered to be dominant, neutron star mat- 
ter can be approximated by pure neutron matter as a 
first step towards a more realistic description of neutron 
stars. Such an approximation was used in the recent 
study [ll| in the model consideration with the effective 
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nuclear forces. It was shown that the behavior of spin 
polarization of neutron matter in the high density region 
in a strong magnetic field crucially depends on whether 
neutron matter develops a spontaneous spin polarization 
(in the absence of a magnetic field) at several times nu- 
clear matter saturation density, or the appearance of a 
spontaneous polarization is not allowed at the relevant 
densities (or delayed to much higher densities) . The first 
case is usual for the Skyrme forces p^l - l23l | , while the sec- 
ond one is characteristic for the realistic nucleon-nucleon 
(NN) interaction [2343^ . In the former case, a ferromag- 
netic transition to a totally spin polarized state occurs 
while in the latter case a ferromagnetic transition is ex- 
cluded at all relevant densities and the spin polarization 
remains quite low even in the high density region. 

The scenario for the evolution of spin polarization at 
high densities in which the spontaneous ferromagnetic 
transition in neutron matter is absent was considered for 
the magnetic fields up to 10^^ G Such an estimate 
for the limiting value of the magnetic field strength in 
the core of a magnetar is usually obtained from the scalar 
virial theorem |33j based on Newtonian gravity. However, 
the density in the core of a magnetar is so large that the 
effects of general relativity might become of importance. 
Then further increase of the core magnetic field is ex- 
pected above 10^^ G [S^. By comparing with the obser- 
vational X-ray data, it was argued that the interior mag- 
netic field strength can be as large as 10^^ G (ssj. Also, 
it was shown in the recent study [36| that in the core of a 
magnetar the magnetic field strength could reach values 
up to 10^° G, if to assume the inhomogeneous distribu- 
tion of the matter density and magnetic field inside a 
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neutron star, or to allow the formation of a quark core 
in the high-density interior of a neutron star (concerning 
the last point, see also Ref. 37|). Under such circum- 
stances, if to admit the interior magnetic fields with the 
strength H > 10^^ G, a different scenario is possible in 
which a field-induced ferromagnetic phase transition of 
neutron spins occurs in the magnetar core. This idea 
was investigated in the recent article [s^, where it was 
shown within the framework of a lowest constrained vari- 
ational approach with the Argonne Vis NN potential that 
a fully spin polarized state in neutron matter could be 
formed in the magnetic field H > 10^^ G. Note, how- 
ever, that, as was pointed out in Refs. [s^ [s^, in such 
ultrastrong magnetic fields the breaking of the 0(3) ro- 
tational symmetry by the magnetic field results in the 
anisotropy of the total pressure, having a smaller value 
parallel than perpendicular to the field direction. The 
possible outcome could be the gravitational collapse of a 
magnetar along the magnetic field, if the magnetic field 
strength is large enough. Thus, exploring the possibil- 
ity of a field-induced ferromagnetic phase transition in 
neutron matter in a strong magnetic field, the effect of 
the pressure anisotropy has to be taken into account be- 
cause this kind of instability could prevent the formation 
of a fully polarized state in neutron matter. This effect 
was not considered in Ref. [ssj . thus, leaving open the 
possibility of the formation of a fully polarized state of 
neutron spins in a strong magnetic field. The degree of 
spin polarization is an important issue for determining 
the neutrino cross sections in the matter, and, hence, it 
is relevant for the adequate description of the neutrino 
transport and thermal evolution of a neutron star . In 
the given study, we provide a fully self-consistent calcu- 
lation of the thermodynamic quantities of spin polarized 
neutron matter at finite temperature taking into account 
the appearance of the pressure anisotropy in a strong 
magnetic field. We consider spin polarization phenom- 
ena in a degenerate magnetized system of strongly inter- 
acting neutrons within the framework of a Fermi liquid 
formalism (40| - |43| , unlike to the previous works [H, HI] , 
where interparticle interactions were disregarded. 

Note that recently new parametrizations of Skyrme 
forces were suggested, BSkl9-BSk21 [11], aimed to avoid 
the spontaneous spin instability of nuclear matter at den- 
sities beyond the nuclear saturation density for the case 
of zero temperature. This is achieved by adding different 
density-dependent terms to the standard Skyrme inter- 
action. The BSkl9 parametrization was constrained to 
reproduce the equation of state (EoS) of nonpolarized 
neutron matter [45| obtained in variational calculation 
with the use of the realistic Urbana Vi4 NN potential and 
the three-body force called there TNI. The BSk20 force 
corresponds to the stiffer EoS [i^, obtained in varia- 
tional calculation with the use of the realistic Argonne 
Vi8 two-body potential and the semiphenomenological 
UIX* three-body force which includes also a relativis- 
tic boost correction. Even a stiffer neutron matter EoS 
was suggested in the Brueckner-Hartree-Fock calculation 



of Ref. [47| based on the same Vis two-body potential 
and a more realistic three-body force containing differ- 
ent meson-exchange contributions. This EoS is the un- 
derlying one for the BSk21 Skyrme interaction. The ad- 
vantage of all of these newly developed Skyrme forces 
is that they preserve the high-quality fits to the mass 
data obtained with the conventional Skyrme forces. An 
important quantity allowing one to distinguish between 
the different representatives of a generalized Skyrme in- 
teraction is the symmetry energy defined as the differ- 
ence between the energies per nucleon in neutron matter 
and symmetric nuclear matter (an alternative definition 
of the symmetry energy is also discussed in Ref. [3). 
In the high density region, the symmetry energy de- 
creases with density for the BSkl9 force, while it in- 
creases with density for BSk20 (moderately) and BSk21 
(steeply) forces. As was clarified in Ref. [48| by testing 
almost 90 parametrizations of the conventional Skyrme 
forces, the Skyrme interactions, predicting the increas- 
ing behavior of the symmetry energy with density, give 
neutron star models in a broad agreement with observa- 
tions (e.g., providing satisfactory description of the mini- 
mum rotation period, gravitational mass-radius relation, 
and the binding energy, released in supernova collapse). 
Considering, based on these arguments, as a more realis- 
tic scenario that in which the symmetry energy increases 
with density in the high density region, in this study we 
will choose the BSk20 Skyrme parametrization for car- 
rying out numerical calculations. Nevertheless, as em- 
phasized in Ref. 4J], only direct experimental evidence 



related to the high densities will allow one to ultimately 
decide which of the BSkl9-BSk21 parametrizations of a 
generalized Skyrme interaction is more appropriate for 
the description of neutron-rich nuclear systems of astro- 
physical interest. 

At this point, it is worthy to note that we consider ther- 
modynamic properties of spin polarized states in neutron 
matter in a strong magnetic field up to the high den- 
sity region relevant for astrophysics. Nevertheless, we 
take into account the nucleon degrees of freedom only, 
although other degrees of freedom, such as pions, hyper- 
ons, kaons. or even quarks could be important at such 
high densities. 



II. BASIC EQUATIONS 

The normal (nonsupcrfluid) states of neutron matter 
are described by the normal distribution function of neu- 
trons /kiK2 = Tr paj^a^j . where k = (p, cr), p is momen- 
tum, a is the projection of spin on the third axis, and 
g is the the density matrix of the system (2l| - [23j . The 
energy of the system is specified as a functional of the 
distribution function /, E = E[f), and determines the 
single particle energy 



£kiK2 (/) — 



dE{f) 
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The self-consistent matrix equation for determining the 
distribution function / foUows from the minimum condi- 
tion of the thermodynamic potential [Z^, and is 



/ = {exp(ro£ -f 
= {exp(yoe) 



-^^.cr, +^4)4-1}" 



(2) 



Here the quantities e, Yi and Y4 are matrices in the space 
of K variables, with (yiA)^^^^ = ^i,4'^«i«2: ^0 = ^/T, 
Yi ~ —Hi/T and Y^ ~ ~iiq/T being the Lagrange mul- 
tipliers, /iQ being the chemical potential of neutrons, and 
T the temperature. In Eq. ©, /x„ = -1.9130427(5)^Ar 
is the neutron magnetic moment (49j (/xjv being the nu- 
clear magneton), (jj are the Pauli matrices. Note that, 
unlike to Refs. [ill, H^, the term with the external mag- 
netic field H is not included in the single partice energy e 
but is separately introduced in the exponent of the Fermi 
distribution 

Further it will be assumed that the third axis is di- 
rected along the external magnetic field H. Given the 
possibility for alignment of neutron spins along or op- 
posite to the magnetic field H, the normal distribution 
function of neutrons and single particle energy e can be 
expanded in the Pauli matrices Ui in spin space 



/(P) = ./o(p)o-0 + /3(p)cr3, 

£(p) = <^o(p)cto + e3(p)CT3. 



(3) 



Using Eqs. © and (U)), one can express evidently the 
distribution functions /o, in terms of the quantities e: 



/o = \{n{i^+)+n{uj-)}, 
h = ^{n{i^+) -n{u^)}. 
Here n{u!) = {exp(ybw) + 1}^^ and 

= £0 — Mo, C3 = —^J-7iH + £3. 



(4) 
(5) 

(6) 



The quantity lu± , being the exponent in the Fermi dis- 
tribution function n, plays the role of the quasiparticle 
spectrum. The branches uj± correspond to neutrons with 
spin up and spin down, respectively. 

The distribution functions / satisfy the normalization 
conditions 
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Here g = g-\- + is the total density of neutron matter, 
g-f and g^ are the neutron number densities with spin 
up and spin down, respectively. The quantity Ag may 
be regarded as the neutron spin order parameter which 
determines the magnetization of the system M = finAg. 



The spin ordering of neutrons can also be characterized 
by the spin polarization parameter 



H 



g 



The magnetization may contribute to the internal mag- 
netic field B = H + 47rM. However, we will assume, 
analogously to the previous studies [MIHlIlSl' that, be- 
cause of the tiny value of the neutron magnetic moment, 
the contribution of the magnetization to the inner mag- 
netic field B remains small for all relevant densities and 
magnetic field strengths, and, hence. 



B « H. 



(9) 



In order to get the self-consistent equations for the 
components of the single particle energy, one has to set 
the energy functional of the system. It represents the 
sum of the matter and field energy contributions 



E{f,H)=E.^{f)+EfiH), Ef{H) 
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The matter energy is the sum of the kinetic and Fermi- 
liquid interaction energy terms [H, [23l | 



E„^if) = Eoif) + E,^tif), 
i?o(/) = 2^£o(p)/o(p). 



(11) 



E^ntif) = Ei=~0(P)/0(P) + £3(P)/3(P)}, 

p 



where 



q 

^~3(p) = :^Et^r(k)/3(q)- 



(12) 
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Here £o(p) = is the free single particle spectrum, 
mo is the bare mass of a neutron, [/,y(k), C/{'(k) are the 
normal Fermi liquid (FL) amplitudes, and £0,63 are the 
FL corrections to the free single particle spectrum. Using 
Eqs. (IlJ and pT|) . we get the self-consistent equations for 
the components of the single particle energy in the form 



Co(p) £o(p) + £o(p) - Aio, 6(P) = -^J■nH 



Hp)- 

(14) 



Taking into account expressions (|4]) and (O for the 
distribution functions /o and /s, solutions of the self- 
consistent Eqs. ((Ti)) should be found jointly with the nor- 
malization conditions (O, ([5]). 

The pressures (longitudinal and transverse with re- 
spect to the direction of the magnetic field) in the system 
are related to the diagonal elements of the stress tensor 
whose explicit expression reads [sij 
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Here 



f-f 



(15) 
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jn = ^{E - TS) - HM is the Helmlioltz free energy 
density, and the entropy S is given by the formula 



5"=-^ ^ {ri(a;cr)lnn(cj„) 
p ""=+, - 
+71(0;^) Inri(ajCT)}, n{uj) = 1 — n(a;). 



(17) 



For the isotropic medium, the stress tensor (jlSp is sym- 
metric. The transverse pt and longitudinal pi pressures 
are determined from the formulas 



Pt 



-CTii = — CT22, Pi = —0-33. 



Hence, using Eqs. ([TU]), (|15p . one can get 



Pt = e 
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where /m = ^ (£",„ — TS") is the matter free energy den- 
sity, and we disregarded in Eqs. ([T5|) . ([T^ the terms pro- 
portional to A/. The structure of the pressures pt and pi 
is different that reflects the breaking of the rotational 
symmetry in the magnetic field. In ultrastrong mag- 
netic fields, the quadratic on the magnetic field term (the 
Maxwell term) will be dominating, leading to increasing 
the transverse pressure and to decreasing the longitudi- 
nal pressure. Hence, at some critical magnetic field, the 
longitudinal pressure will vanish, resulting in the longi- 
tudinal instability of neutron matter. Obviously, at finite 
temperature the pressures pt and pi will be larger com- 
pared to the zero temperature case, and, hence, increase 
of the temperature will lead to the increase of the critical 
magnetic field. Here we would like to find the magnitude 
of the critical field at temperatures of about a few tens of 
MeV, which can be relevant for protoneutron stars, and 
also to determine the corresponding maximum degree of 
spin polarization in neutron matter. 



temperature a spontaneous ferromagnetic phase transi- 
tion could occur. Actually, this is the case as will be 
shown later. In the model calculations of this section 
we consider the temperatures somewhat larger than the 
temperatures which could be reachable in the interior of 
protoneutron stars [Hj. This will help us to find the 
critical temperature above which a spontaneous polar- 
ization appears, and will also allow us to determine the 
relevant temperature range for studying spin polarization 
at iJ 0. 

The recently developed parametrizations BSkl9- 
BSk21 of the Skyrme effective forces appear as a gener- 
alization of Skyrme effective NN interaction of the con- 
ventional form. In the conventional case, the amplitude 
of Skyrme NN interaction reads j53| 



v{p,q) = toil + XoPa) + Itsil + X3P^)g'' 
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(20) 



— ti(l + xiF,)(p2 



— (1 +X2P^)pq, 



where = {1 + (7\<Ji)l2 is the spin exchange operator, 
tj, Xi and a. are some phenomenological parameters spec- 
ifying a given parametrization of the Skyrme interaction. 
The Skyrme interaction used in Ref. [4J| has the form 



^^'(P,q) =i'(p,q) + 1^*4(1 + a;4P.)(p' + q') (21) 



ft2 



^5(1 + a;5/'cr)pq. 



In Eq. (PT|) . two additional terms are the density- 
dependent generalizations of the t\ and ti terms of the 
usual form. Specific values of the parameters ti,Xi,a,P 
and 7 for Skyrme forces BSkl9-BSk21 are given in Ta- 
ble 1 [3. 

The normal FL amplitudes Uo,Ui can be expressed in 
terms of the Skyrme force parameters. For conventional 
Skyrme force parametrizations, their explicit expressions 
are given in Refs. [4l|,|4^. As follows from Eqs. (|20p and 
(|2ip , in order to obtain the corresponding expressions for 
the generalized Skyrme interaction (PT|) . one should use 
the substitutions 

ti ^ ti+ t^g^^, tixi ^ tixi + t4X4g^, (22) 
^2^*2+ he'', t2X2 t2X2 + hx^g''. (23) 

Therefore, the FL amplitudes are related to the pa- 
rameters of the Skyrme interaction (|2ip by formulas |54{ 



III. SPIN POLARIZATION AT // = 0, T / 

For providing numerical calculations, we use the BSk20 
Skyrme interaction (43 | developed to reproduce the zero 
temperature microscopic EoS of nonpolarized neutron 
matter (46j . Although spontaneous spin polarization at 
zero temperature is missing for this parametrization for 
all relevant densities, it is not excluded that at finite 



(7o"(k) =2to(l-.To) + |e"(l-a;3) 
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(24) 

+ ^4(1 - X4)g^ + 3t2(l + X2) + 3^5(1 + X5)g'^]k'^, 

Ul\k) = -2io(l - xo) - |e"(l - .T3) + ^[t2{l + X2) 

(25) 

+ h{l + X5)^)^ - ti(l - a-i) - U{1 - X4)/]k2 
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TABLE I. The parameters of the BSkl9-BSk21 Skyrme forces, 
according to Ref. [13]. The value of the nuclear saturation 
density go is shown in the bottom line. 







BSkig 


BSk20 


BSk21 


to 


[MeV fm^] 


-4115.21 


-4056.04 


-3961.39 


tl 


[MeV fm^] 


403.072 


438.219 


396.131 




[MeV fm^] 











ts 


[MeV fm='+^"] 


23670.4 


23256.6 


22588.2 


U 


[MeV fm^+^'^] 


-60.0 


-100.000 


-100.000 


tz, 


[MeV fm^+^^l 


-90.0 


-120.000 


-150.000 


Xo 




0.398848 


0.569613 


0.885231 


Xl 




-0.137960 


-0.392047 


0.0648452 


t2X2 [MeV fm^] 


-1055.55 


-1147.64 


-1390.38 


X3 




0.375201 


0.614276 


1.03928 


X4 




-6.0 


-3.00000 


2.00000 


X5 




-13.0 


-11.0000 


-11.0000 


a 




1/12 


1/12 


1/12 


P 




1/3 


1/6 


1/2 


7 




1/12 


1/12 


1/12 


Qo 


[l/fm^] 


0.1596 


0.1596 


0.1582 



Now wc present the results of the numerical solution 
of the self-consistent equations at if = with the BSk20 
Skyrme force. Fig. [1] shows the spin polarization param- 
eter of neutron matter as a function of the density for 
a few fixed values of the temperature of about several 
tens of MeV. At zero temperature, there is no sponta- 
neous polarization at all relevant densities because two 
additional terms in a generalized form ([21]) of the Skyrme 
interaction were constrained just with the aim to exclude 
a nonzero polarization at vanishing temperature. Spon- 
taneous polarization docs not appear up to some critical 
temperature which is, at least, larger than 35 MeV. 
Beyond T^, spontaneous spin polarization exists in a fi- 
nite density interval [qci, Qc^)- The unexpected moment 
is that the temperature promotes spontaneous spin polar- 
ization increasing both the width of the density domain 
where a nonzero polarization exists and the magnitude 
of the spin polarization parameter. In particular, if to 
approach the density interval {QcnQc2) from the lower 
densities then the left critical point , at which sponta- 
neous polarization appears, decreases with temperature, 
contrary to intuition which suggests that the tempera- 
ture should act as a preventing factor to spin polarization 
and, hence, should delay its appearance. Analogously, 
with increasing temperature the right critical point 
for the disappearance of spontaneous polarization should, 
according to intuition, decrease, contrary to what really 
occurs, i.e., the critical density increases with tem- 
perature. 

In order to clarify whether a spontaneously spin polar- 
ized state is thermodynamically preferable over the non- 
polarized state, one should compare the corresponding 
free energies. Fig.[2]shows the difference between the free 
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FIG. 1. Neutron spin polarization parameter as a function of 
the density for the BSk20 Skyrme force at -ff = and several 
fixed values of the temperature. 



energies per neutron of spin polarized and nonpolarized 
states, 5F/A = {F{g,T,n{g,T)) - F{g,T,n = 0))/A, 
as a function of the density at the same fixed tempera- 
tures considered above. It is seen that a spontaneously 
polarized state is preferable over the nonpolarized state 
for all relevant densities and temperatures where spon- 
taneous polarization exists. With increasing the tem- 
perature, the minimum of the difference 5F / A becomes 
more pronounced, and, hence, a spontaneously polarized 
state becomes more stable with respect to the nonpolar- 
ized one. Thus, the state with spontaneous polarization, 
described by the spin polarization parameter with such 
unusual properties (cf. Fig. [T]), is supported thermody- 
namically by the balance of the free energies. 

In order to get a deeper insight into the problem, let us 
consider separate contributions to the difference between 
the free energies per neutron 6F/A = SE/A — T5S/A. 
Fig. [3] shows the difference between the energies per neu- 
tron of spin polarized and nonpolarized states, SE/A = 
{E{Q,T,U{g,T)) - E{Q,T,n = 0))/A, as a function of 
the density at the same fixed temperatures considered 
above. It is seen that the energy per neutron of a spin 
polarized state is always larger than that of the nonpolar- 
ized state for the density domain where spontaneous spin 
polarization exists. This is because increasing the tem- 
perature and spin polarization leads to increasing the ki- 
netic energy term in the energy functional of the system. 
The sign of the difference SE/A could be, in principle, 
inverted by the negative contribution of the term in the 
energy functional describing spin correlations in neu- 
tron matter with nonzero polarization, but that is not, 
however, the case. Therefore, the inequality SF/A < 
can hold only because of the inequality SS/A > for 
the density range where spontaneous polarization exists. 
Fig. m shows that this is actually true, and the entropy 
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FIG. 2. The difference between the free energies per neutron 
of spin polarized and nonpolarized states as a function of the 
density at = and several fixed values of the temperature 
for the BSk20 Skyrme force. The difference is shown only for 
the density domains where spontaneous polarization exists. 



per neutron of a spin polarized state is larger than that 
of the nonpolarized state for the corresponding temper- 
atures and densities. This unexpected behavior is con- 
tradicting to intuition which suggests that the entropy of 
a more ordered spin polarized state should be less than 
that of the nonpolarized state. Note that such an un- 
usual behavior of the entropy of a spin polarized state 
was found earlier for neutron matter with the Skyrme ef- 
fective interaction (sBj and for symmetric nuclear matter 
with the Gogny effective interaction [sl, [s^] (in the lat- 
ter case, for antiferomagnetically ordered nucleon spins) . 
The difference, however, is that in these earlier studies 
instability with respect to spontaneous spin ordering oc- 
curred already at zero temperature whereas in the given 
case it appears only at temperatures lar ger than the crit- 
ical one. Also, it was clarified earlier jSSl [57l | that the 
unusual behavior of the entropy of a spin polarized state 
should be traced back to its dependence on the effective 
masses of spin-up and spin-down nucleons and to a viola- 
tion of a certain constraint on them at the corresponding 
temperatures and densities. In Ref. [s^, this constraint 
was formulated for a totally polarized neutron matter, 
and in Ref. 



57| for symmetric nuclear matter with arbi- 



trary antiferromagnetic spin polarization. 

Let us verify now whether this holds true in our case. 
In the low-temperature limit the entropy per neutron is 
given by expression 



2£Fa 



-T, 



(26) 



2 , 2 

FtT 



h k • 

where £„ = ^m" is the Fermi energy of neutrons with 

spin up and spin down, and fco- = (GTr^pg-)^^"^ is the re- 
spective Fermi momentum. The low-temperature expan- 



FIG. 3. Same as in Fig. [2] but for the difference between 
the energies per neutron of spin polarized and nonpolarized 
states. 



sion (|26p is valid till T/epa ^ 1- By requiring for the 
difference between the entropies of spin polarized and 
nonpolarized states to be negative, one can derive the fol- 
lowing constraint on the effective masses m„-|- and m„4, of 
neutrons with spin up and spin down in a spin polarized 
state IH: 



Writ , 

D = ^(1 + n) 



(l-n)3^2<0, (27) 



where 



2m^(;) 2m 



+ ^-^[t2{l + X2)+h{l + x,)e''] (28) 

^ 



+ ^[ii(l-xi)+t4(l-a;4)/ 

In the constraint (P7)) . the effective mass m„ of a neu- 
tron in nonpolarized neutron matter is given by (s^ 



+ |[<i(l-xi)+t4(l-a;4)/ (29) 



2?Ti„ 2mo 8 

+ it2{l+X2)+2,t5{l+Xz)Q^]. 

After the self-consistent determination of the spin po- 
larization parameter, one can check whether the inequal- 
ity p7)) is satisfied at the corresponding densities and 
temperatures. Fig. [5] shows the left-hand side D of the 
constraint ([27|) for the branch n(gi, T) of spontaneous po- 
larization as a function of the density at the temperatures 
T = 37 MeV and T = 40 MeV, at which the accuracy of 
the approximation T/epa ^ 1 is satisfactory. It is seen 
that the inequality (P7)) is violated implying that the en- 
tropy of a spontaneously polarized state is larger than the 



7 



0.40 



0.35 - 



o 0.30 



CO 0.15 

to 



0.05 
0.00 







T = 37 MeV 


- 


/' 


40 MeV 




/ 


^\ 60 MeV 


- 


/ 


\ 




/ 


\ 




/ 

/ 

/ 


\ 

\ 

\ 


- 


; 

/ 


\ 

\ 


/ 




^- BSk20 


/ 

/ 


/ 


\ H = 


( 


/ 

/ 






/ 

/ 




il 




■ — ^ ^ 



1.0 1.5 2.0 2.5 3.0 3.5 

p/po 

FIG. 4. Same as in Fig. [2] but for the difference between 
the entropies per neutron of spin polarized and nonpolarized 
states. 



entropy of the nonpolarized state at the respective den- 
sities and temperatures. Hence, the unusual behavior of 
the entropy of a spontaneously polarized state mentioned 
above can be related to the peculiarities of its dependence 
on the effective masses of neutrons with spin up and spin 
down. A nontrivial character of the density dependence 
of the effective masses m„-f and m^j, in neutron matter 
with spontaneous polarization at different temperatures 
is clearly seen from Fig. [6] 

In the subsequent analysis, following the scenario ac- 
cording to which spontaneous polarization should be 
avoided at the relevant densities and temperatures, we 
will confine our analysis to the temperatures up to 30 
MeV which are, definitely, less than the critical tempera- 
ture Tc > 35 MeV. Such a choice of the relevant tempera- 
ture interval is consistent with the results of a completely 
independent research (52j of hybrid stars in the context 
of relativistic mean-field theory, according to which the 
maximum temperature attainable in their interior does 
not exceed 35 MeV. 



0.25 
0.20 
0.15 
0.10 
0.05 
0.00 



T = 37 MeV 
40 MeV 



BSk20 
H = 




1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 

p/po 
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37 MeV and T = 40 MeV for the BSk20 Skyrme force. 
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IV. LONGITUDINAL AND TRANSVERSE 
PRESSURES AT FINITE TEMPERATURE. 
ANISOTROPIC EOS 

In this section, we will study the influence of finite tem- 
peratures on thermodynamic quantities of spin polarized 
neutron matter in an ultrastrong magnetic field. We will 
take into account the effects of the pressure anisotropy, 
and. in particular, will clarify to which extent the crit- 
ical magnetic field, at which the longitudinal instability 
in magnetized neutron matter occurs, will increase due 
to the impact of finite temperatures. 

First, we present the results of the numerical solution 



of the self-consistent equations. Fig. [7] shows the spin po- 
larization parameter of neutron matter as a function of 
the magnetic field H at two different temperatures, T = 
and T = 30 MeV, and at two different values of the neu- 
tron matter density, q = Sgo and q = 4po, which can 
be relevant for the central regions of a magnetar. Under 
increasing the density, the effect produced by the mag- 
netic field on spin polarization of neutron matter becomes 
smaller. It is seen that the impact of the magnetic field 
remains insignificant up to the field strength H ~ 10^'' G. 
At the magnetic field H = 10^® G, usually considered 
as the maximum magnetic field strength in the core a 
magnetar (according to a scalar virial theorem [33|). the 



8 





" "^"^^"^^-^^^^ T = OMeV 




" r\ \ 30 MeV 




^ \ \ 
\ \ \ 

^ \ \ 
\ \ \| 

\ \ 




\\\^P = 4Po 




^\\ 

P = 3P„ \\ \ 

\ \ \ 




\ \ \ 

A \ BSk20 

W. - 



oi — — ^- V . I 

10'' 10'" 10'' 10'" 10' 

H[G] 



FIG. 7. Neutron spin polarization parameter as a function 
of the magnetic field strength for the BSk20 Skyrme force at 
T = and T — 30 MeV, and at two fixed densities, g = 3qo 
and g = Ago. The vertical arrows indicate the maximum 
magnitude of spin polarization attainable at the given tem- 
perature and density, see further details in the text. 

magnitude of the spin polarization parameter doesn't ex- 
ceed 45% a.t g — 3po and 19% a.t g ~ igo (for the tem- 
peratures under consideration). However, the situation 
changes if the larger magnetic fields are allowable: With 
further increasing the magnetic field strength, the mag- 
nitude of the spin polarization parameter increases and 
spin polarization approaches its limiting value 11 = —1, 
corresponding to a fully spin polarized state. For exam- 
ple, a fully polarized state is formed at _ff w 1.3 • 10^^ G 
for the temperature T = MeV and at H ^ 2.3 • 10^^ G 
for T = 30 MeV at g = 3 go, i.e., certainly, for mag- 
netic fields H > 10^^ G. Note that we speak about a 
fully polarized state at finite temperature although some 
quantity of neutrons with spin up are always present at 
T ^ 0. Nevertheless, this quantity may be made arbi- 
trary small with further increasing the magnetic field, 
and we consider that a fully polarized state is formed, 
if the deviation from the limiting value 11 = —1 is less 
than 10~^. With increasing the temperature, the value 
of the magnetic field, at which a fully polarized state oc- 
curs, increases, as one could expect. However, practically 
up to magnetic fields of about 10^^ G, spin polarization 
demonstrates the unusual behavior and increases with 
temperature. Further it will be shown that this behav- 
ior is thermodynamically supported by the corresponding 
balance of the Helmholtz free energies. The meaning of 
the vertical arrows in Fig. [7] is explained later in the text. 

Now, we should check whether a fully spin polarized 
state of neutrons in a strong magnetic field can indeed be 
formed by calculating the anisotropic pressure in dense 
neutron matter. Fig. [8^ shows the pressures (longitu- 
dinal and transverse) in neutron matter as functions of 
the magnetic field H at the same fixed temperatures and 



densities, considered above. The upper branches in the 
branching curves correspond to the transverse pressure, 
the lower ones to the longitudinal pressure. First, it is 
clearly seen that up to some threshold magnetic field 
the difference between the transverse and longitudinal 
pressures is unessential that corresponds to the isotropic 
regime. Beyond this threshold magnetic field strength, 
the anisotropic regime holds for which the transverse 
pressure increases with H while the longitudinal pres- 
sure decreases. The increase of the temperature leads 
to the increase of the pressures, transverse pt and lon- 
gitudinal pi. Also, the increase of the density has the 
same effect on the pressures pt and pi as the increase of 
the temperature. The most important feature is that the 
longitudinal pressure vanishes at some critical magnetic 
field He marking the onset of the longitudinal instability 
in neutron matter. For example. He ~ 1.56 • 10^* G for 
T = MeV and He « 1.64 • 10^® G for T 30 MeV at 
g = 3(?o, and H^ « 2.42 ■ 10^** G for T MeV and 
He ~ 2.48 ■ 10^*^ G for T 30 MeV at (? = 4,go. Hence, 
at finite temperatures relevant for protoneutron stars the 
critical magnetic field is increased compared to the zero 
temperature case but this increase is, in fact, insignifi- 
cant. Even with accounting for the finite temperature 
effects, the critical field doesn't exceed 10^^ G for the 
density range under consideration. 

The magnitude of the spin polarization parameter H 
cannot also exceed some limiting value corresponding to 
the critical field He- These maximum values of the H's 
magnitude are shown in Fig. [7] by the vertical arrows. In 
particular, H^ « -0.46 for T = MeV and H^ « -0.58 
for T = 30 MeV at g = 3go, and H^ « -0.38 for 
T = MeV and H^ « -0.41 for T = 30 MeV at p = Ago. 
As can be inferred from these values, the appearance of 
the negative longitudinal pressure in an ultrastrong mag- 
netic field prevents the formation of a fully polarized spin 
state in the core of a magnetar. Therefore, only the on- 
set of a field-induced ferromagnetic phase transition, or 
its near vicinity, can be caught under increasing the mag- 
netic field strength in dense neutron matter at finite tem- 
perature. A complete spin polarization in the magnetar 
core is not allowed by the appearance of the negative 
pressure along the direction of the magnetic field, con- 
trary to the conclusion of Ref. [s^ where the pressure 
anisotropy in a strong magnetic field was disregarded. 

Fig. [HJd shows the difference between the transverse 
and longitudinal pressures normalized to the value of the 
pressure po in the isotropic regime (which corresponds to 
the weak field limit with pi = pt = Po)'- 

Po 

Applying for the transition from the isotropic regime to 
the anisotropic one the criterion 5 ~ 1. the transition 
occurs at the threshold field Hth ~ 1.15 ■ 10^^ G for 
T = MeV and Hth ~ 1-22 • 10^*^ G for T = 30 MeV 
at g = 3go, and at Hth ~ 1-83 • lO^^ G for T = MeV 
and Hth « 1.86 • lO^^ G for T = 30 MeV at g ^ igo- 
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neutron matter. 



In all cases under eonsideration, the threshold field Hth 
is greater than 10^* G, and, hence, the isotropic regime 
holds for the fields up to 10^^ G. For comparison, the 
threshold field for a relativistic dense gas of free charged 
fermions at zero temperature was found to be about 
10^^ G [11] (without including the anomalous magnetic 
moments of fermions). For a degenerate gas of free neu- 
trons at zero temperature the model dependent estimate 
gives Hth ^ 4.5 ■ 10^* G [H] (including the neutron 
anomalous magnetic moment). The normalized split- 
ting of the transverse and longitudinal pressures increases 
more rapidly with the magnetic field at the smaller den- 
sity and/or at the lower temperature. The vertical ar- 
rows in Fig. HId indicate the points corresponding to the 
onset of the longitudinal instability in neutron matter. 
Since the threshold field Hth is less than the critical field 
He for the appearance of the longitudinal instability, the 
anisotropic regime can be relevant for the core of a mag- 
netar. The maximum allowable normalized splitting of 
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FIG. 9. (a) Same as in Fig. [7] but for the Helmholtz free 
energy density of the system, (b) Same as in Fig. [7] but for 
the ratio of the magnetic field energy density to the Helmholtz 
free energy density of the system. The meaning of the vertical 
arrows is the same as in Fig. [SJd. 



the pressures corresponding to the critical field He is 
6 2. If the anisotropic regime sets in, a neutron star 
has the oblate form. Thus, as follows from the preced- 
ing discussions, in the anisotropic regime the pressure 
anisotropy plays an important role in determining the 
spin structure and configuration of a neutron star. 

At the given thermodynamic variables g, T and H, the 
Helmholtz free energy is a relevant thermodynamic func- 
tion, whose minimum determines a state of thermody- 
namic equilibrium. Fig. shows the Helmholtz free 
energy density of the system as a function of the mag- 
netic field H at two fixed temperatures, T — and 
T = 30 MeV, and at two different densities, g = 3go 
and g = Ago- It is seen that the magnetic fields up 
to H ^ 10^® G have practically small effect on the 
Helmholtz free energy density fn, but beyond this field 
strength the contribution of the magnetic field energy to 
the free energy fn rapidly increases with H. However, 
this increase is limited by the values of the critical mag- 
netic field corresponding to the onset of the longitudinal 
instability in neutron matter. The respective points on 
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the curves are indicated by the vertical arrows. 

Fig. |9b shows the ratio of the magnetic field energy 
density ef — ^ to the Helmholtz free energy density 
under the same assumptions as in Fig. [5^. The in- 
tersection points of the respective curves in this panel 
with the line ef/fu = 0.5 correspond to the magnetic 
fields at which the matter and field contributions to the 
Helmholtz free energy density arc equal. This happens 
at iJ « l.lS-lQis GforT = OMeV and i? w LOS-IO^^ G 
for T = 30 McV &i g = Spo, and at « 1.81 • IQ^^ G for 
T = MeV and H w 1.76 • 10^^ G for T = 30 MeV at 
Q = 4£io- These values are quite close to the respective 
values of the threshold field Hth, and, hence, the tran- 
sition to the anisotropic regime occurs at the magnetic 
field strength at which the field and matter contributions 
to the Helmholtz free energy density become equally im- 
portant. It is also seen from Fig.[3]3 that in all cases when 
the longitudinal instability occurs in the magnetic field 
He the contribution of the magnetic field energy density 
to the Helmholtz free energy density of the system dom- 
inates over the matter contribution. 

Because of the pressure anisotropy, the EoS of neu- 
tron matter in a strong magnetic field is also anisotropic. 
Fig. [TU] shows the dependence of the Helmholtz free en- 
ergy density fn on the transverse pressure (top panel) 
and on the longitudinal pressure (bottom panel) after ex- 
cluding the dependence on H in these quantities. Since 
the dominant Maxwell term enters the pressure pt and 
free energy density Jh with positive sign and the pres- 
sure pi with negative sign, the free energy density fu 
is the increasing function of pt and decreasing function 
of pi. In the case of fnipt) dependence, at the given 
density, the same pt corresponds to the larger magnetic 
field H at the temperature T — MeV compared to the 
T = 30 MeV case (sec Fig. Hi). The overall effect of 
two factors (temperature and magnetic field) will be the 
larger value of the free energy density fn at the given pt 
and density for the temperature T = MeV compared 
with the T = 30 MeV case (see Fig. The analo- 

gous arguments show that, at the given temperature and 
Pt, the Helmholtz free energy density is larger for the 
smaller density. In the case of fnipi) dependence, at 
the given density, the same pi corresponds to the smaller 
magnetic field H for the temperature T = MeV com- 
pared to the T = 30 MeV case (see Fig. [5^). Hence, 
the free energy density fn at the given pi and density 
is larger for the temperature T = 30 MeV than that for 
the T = MeV case (see Fig. ^Ujp). Analogously, at the 
given temperature and pi, the free energy density fn is 
larger for the larger density. In the bottom panel, the 
physical region corresponds to the positive values of the 
longitudinal pressure. 

It is worthy to note at this point that since the EoS 
of neutron matter becomes essentially anisotropic in an 
ultrastrong magnetic field, the usual scheme for find- 
ing the mass-radius relationship based on the Tolman- 
Oppenheimer-Volkoff (TOV) equations [HI for a spher- 
ically symmetric and static neutron star, should be re- 
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FIG. 10. The Helmholtz free energy density of the system as 
a function of: (a) the transverse pressure pt, (b) the longitudi- 
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lines), and at two fixed densities, q = Sfio and q — Aqq. The 
meaning of the vertical arrows in the top panel is the same 
as in Fig.[8}D. In the bottom panel, the physical region corre- 
sponds to pi > 0. 



vised. Instead of this, the corresponding relationship 
should be found by the self-consistent treatment of the 
anisotropic EoS and axisymmetric TOV equations sub- 
stituting the conventional TOV equations in the case of 
an axisymmetric neutron star. 



V. UNUSUAL BEHAVIOR OF THE ENTROPY 

AT H ^0 

As was discussed in the previous section, the magni- 
tude of the spin polarization parameter increases with 
temperature in the fields up to about 10^^ G. The 
Helmholtz free energy density fn , whose minimum at the 
given g, T, H determines the state of a thermodynamic 
equilibrium, decreases with temperature (cf. Fig. [9^) 
and. hence, such an unusual behavior of spin polariza- 
tion with temperature is supported thermodynamically. 
The Helmholtz free energy density fn can be decomposed 
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into the matter and field contributions, 

fn = fHm + 6/ , 

with the matter contribution being fHm — ^(-Em— T5) — 
HM. The decrease of the Helmholtz free energy with 
temperature is, therefore, to be attributed to its matter 
part. Fig. [TT] explicitly shows this point. 

An unexpected moment appears if we consider sep- 
arately the behavior of the entropy of neutron matter 
with a generalized Skyrme interaction in a strong mag- 
netic field. In Fig.[T21 the difference between the entropy 
per neutron of magnetized neutron matter and that of 
the nonpolarized state (with 11 = at if = 0) is pre- 
sented as a function of magnetic field at the tempera- 
tures T = 15 MeV and T = 30 MeV, and at the same 



densities regarded above. It is seen that this difference 
is positive for all relevant magnetic field strengths. It 
looks like a spin polarized state is less ordered than the 
nonpolarized one, contrary to intuitive assumption. In 
section IIIIl we showed that the unusual behavior of the 
entropy of a spontaneously polarized state is related to 
its dependence on the effective masses of neutrons with 
spin up and spin down and to the violation of the cri- 
terion pT]). The entropy of magnetized neutron matter 
is given by the same general expression (jl7|) . and, after 
providing the low-temperature expansion, we would ar- 
rive at the same constraint (|27p on the effective masses 
in a spin polarized state guaranteeing that its entropy is 
less than that of the nonpolarized state. Fig. [T3] shows 
the left side D of the constraint (|27p as a function of the 
magnetic field strength at the temperature T = 15 MeV, 
and densities g = 3go and g = Ago, at which the accuracy 
of the approximation T/epa ^ 1 is acceptable. It is seen 
that the criterion ^7} is violated, and this explains the 
unusual behavior of the entropy of dense neutron matter 
in a strong magnetic field shown in Fig. 1121 

Note that the unconventional behavior of the entropy 
of magnetized neutron matter with Skyrme interaction 
was found earlier in Ref. [HO]. The difference is that for 
the SLyT Skyrme interaction used in that work a spon- 
taneously polarized state appears already at zero tem- 
perature, while in the given research with a newly devel- 
oped BSk20 Skyrme force spontaneous polarization ap- 
pears only at temperatures above the critical one. We 
have checked that the last feature is also characteristic 
for the BSkl9 and BSk21 Skyrme forces. If to consider 
the appearance of a spontaneously polarized state as a 
weak point of a certain Skyrme parametrization (just 
this argument was used in Ref. [4J] as the motivation 
for developing a new series of Skyrme forces) then this 
underlines the necessity to further concentrate the ef- 
forts on building a new generation of Skyrme forces be- 
ing free of such kind of spin instabilities. Such an at- 
tempt was made in the recent article [59| by attracting 
the ideas from the nuclear energy density functional the- 
ory. However, the constraints obtained in this study on 
the Skyrme force parameters lead to the unrealistic con- 
sequence that the effective masses of nucleons with spin 
up and spin down in a polarized state should be equal, 
contrary to the results of calculations with realistic NN 
interaction [2^ |3l|. On the other hand, the observa- 
tional data still do not rule out the existence of a fer- 
romagnetic hadronic core inside a neutron star caused 
by spontaneous ordering of hadron spins (in this respect, 
see, e.g., Refs. I63,[6l|). In any case, developed recently 
generalized Skyrme parametrizations BSkl9-BSk21 are, 
currently, among the most competitive Skyrme forces for 
providing neutron star calculations, and, certainly, they 
are suitable for getting a qualitative estimate of the ef- 
fects of the pressure anisotropy in strongly magnetized 
neutron matter at finite temperature. 

In summary, we have considered spin polarized states 
in dense neutron matter in the model with the recently 
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developed BSk20 Skyrme interaction at finite tempera- 
ture under the presence of strong magnetic fields up to 
10^° G. Although the BSk20 Skyrme force was worked 
up with the aim to avoid spontaneous spin instability at 
zero temperature, it has been shown that spontaneous in- 
stability appears at temperatures above the critical one, 
which is, at least, larger than 35 MeV. By this reason, 
we limited our consideration by the temperatures up to 
30 MeV. For a spontaneously polarized state at finite 
temperature, the entropy demonstrates the unusual be- 
havior being larger than that of the nonpolarized state. 
This feature has been related to the dependence of the 
entropy of a spin polarized state on the effective masses 
of spin-up and spin-down neutrons and to the violation of 
some constraint on them at the corresponding densities 



and temperatures. In strong magnetic fields considered 
in this study the total pressure in neutron matter be- 
comes anisotropic. It has been shown that for the mag- 
netic fields H > Hth ^ 10^* G the pressure anisotropy 
has a significant impact on thermodynamic properties of 
neutron matter. In particular, vanishing of the pressure 
along the direction of the magnetic field in the critical 
field He > Hth leads to the appearance of the longitu- 
dinal instability of neutron matter. With increasing the 
density and temperature of neutron matter, the thresh- 
old Hth and critical H^ magnetic fields also increase. In 
the limiting case considered in this study and correspond- 
ing to the density of about four times nuclear saturation 
density and the temperature of about a few tens of MeV, 
the critical field He doesn't exceed 10^^ G. This value can 
be considered as the upper bound on the magnetic field 
strength inside a magnetar. Our calculations show that 
the appearance of the longitudinal instability prevents 
the formation of a fully spin polarized state in neutron 
matter, and only the states with moderate spin polar- 
ization can be developed. In the anisotropic regime, the 
field contribution to the Hclmholtz free energy density 
becomes comparable and even dominates over the mat- 
ter contribution. The longitudinal and transverse pres- 
sures and anisotropic EoS of neutron matter in a strong 
magnetic field have been determined at the densities and 
temperatures relevant for the interior of a magnetar. It 
has been clarified that the entropy of strongly magnetized 
neutron matter with the Skyrme BSk20 force demon- 
strates the unusual behavior similar to that of the en- 
tropy of spontaneously polarized state. In both cases, 
the same reason, discussed above, is responsible for such 
a behavior. The obtained results can be of importance 
in the studies of cooling history and structure of strongly 
magnetized neutron stars. 
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